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I.  Introduction 

The  electromagaetic  field  which  can  he  constructed  hy  the 

hi 
methods  of  geometrical  optics,  following  R.Z,  Lunehxirg'-  -^  becomes 

infinite  at  caustics  of  the  optical  ray  systems.  In  this  paper  we 

overcome  this  difficulty  and  obtain  a  finite  value  for  the  field  at 

a  caustic  which  is  produced  by  the  reflection  of  any  incident  wave 

from  an  arbitrary  reflector  (all  in  two  dimensions  for  simplicity). 

What  we  actially  determine  is  the  leading  term  In  the  asymptotic 

expansion  of  the  field  with  respect  to  k  =  r—  ,  for  large  k.  Off 

the  caustic  this  is  just  the  geometrical  optics  field,  but  on  the 

1  _  1 
caustic  it  is  proportional  to  k?   n  where  n  is  an  integer  greater 

than  two.  At  a  perfect  focus,  such  as  the  focus  of  a  parabolic 

1 
cylinder,  the  field  is  proportioiial  to  k^  .  Thus  on  a  caustic  or  at 

a  focus  the  field  increases  indefinitely  as  k  does. 

We  have  exemplified  these  general  results  by  considering 
the  special  cases  of  a  plane  wave  incident  on  the  concave  sides  of 
a  segment  of  a  parabolic  cylinder  and  a  segment  of  a  eirc-ulsT  cylinder. 
In  the  parabolic  case,  with  incidence  along  the  axis  of  the  parabola, 
the  entire  caustic  is  a  point,  the  focus  of  the  parabola.  The 
asymptotic  behavior  of  the  field  at  the  focus  is  obtained  explicitly 
and  transition  formulae  which  show  how  the  field  varies  in  the 
neighborhood  of  the  focus  are  also  obtained.  These  formulae  in- 
volve integrals  which  we  have  evaluated  numerically  in  some  special 
cases. 

In  the  case  of  the  circular  cylinder,  the  caus  iic  is  a 
curve  with  a  cusp.  The  asymptotic  behavior  at  the  cusp  is  different 
from  that  on  the  rest  of  the  caustic.  The  variation  of  the  field 
near  the  caustic  is  described  in  terms  of  Hardy's  generalized  Airy 
functions. 

Since  our  results  yield  exactly  the  leading  term  in  the 
asymritotic  expansion  of  the  field,  they  constitute  a  check  on  the 
solutions  proposed  by  Debye'-  JPichtL  J  and  Luneburg  -^  .  Each 
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of  these  authors  constinicted  e  solution  of  the  wave  equation,  or  of 
IfejDvell's  equations,  which  is  re^rular  i:i  en  infinite  homogeneous  space 
end  satisfies  coniitlons  at  infinity,  Luneturg  proved  that  the  solu- 
tion thus  ohteined  reduces  to  the  kno'Am  geometrical  optics  solution  as 
k  hecomes  infinite.   In  an  actual  prohlem  involving  houndaries  the  cor- 
rect solution  nust  satisfy  certain  'boundary"  conditions.  Since  the  solu- 
tions Just  nientioned  do  not  satisfy  such  conditions,  it  is  not  certain 
that  they  will  yield  the  correct  field  at  a  caustic  in  a  reel  prohlem 
involving  houndcries.  However,  we  have  proved  that  they  do  yield  erect- 
ly  the  sariie  asymptotic  field  on  a  caustic  as  is  given  hy  our  method,  and 
therefore  their  use  in  practical  problems  is  Justified. 
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II.  Fora-glation;  The  Geometrical  Optics  Field 

We  consider  the  solution  E(x,y)  of  the  reduced  wave  equation 

(1)  (7^  +  k^)  E  »  0. 

The  solution  'nu.st  satisfy  the  houndsry  conditions 

(2)  B  =  0   on  C  , 

Here  C  is  a  given  piecewise  smooth  curve  which  we  will  call  a  reflector. 
At  any  corners  of  C  the  derivatives  of  E  may  "become  singular.  Therefore 
we  must  impose  a  regularity  condition  at  these  comers,  ©.s»  that  the 
singularity  he  the  weakest  possihle  one,  or  that  the  "  energy"  integral  "be 
finite.   In  addition  ve  assume  that  an  incident  field  E^(x,y)  which 
satisfies  (l)  is  given,  and  that  the  reflected  field  E  =  E  -  E.  satisfies 
the  radiation  condition 

(3)  E^  =  E  -  E^  is  outgoing. 

These  conditions  complete  the  formulation  of  the  pro"blem.  This  formulation 
applies  to  the  z  coaroonent  of  the  electric  field  when  the  reflector  is  a 
perfect  conductor  and  the  surrounding  medium  is  homogeneous. 

Rather  than  find  E  itself,  we  treat  the  simpler  prohlem  of 
determining  the  leading  term  in  the  asymptotic  expansion  of  E  with 
respect  to  k,  for  large  k.  This  is  what  we  call  the  geometrical  optics 
field.  We  "begin  with  the  leading  term  in  the  asya:ptotic  expansion  of 
the  incident  field  B.,  We  assume  it  to  be  of  the  form 

(^)  \  ^  A^(x,y)  exp  [ib»^  (x,y3  . 

Now,  following  Lunehurg,  we  know  that  the  zero- order  term  in 
the  asymptotic  expansion  of  the  field  is  zero  in  the  shadow.  Further- 
more, the  leading  term  in  the  reflected  field  is 

(5)  \^  A^(x,7)  exp  [ibi;^  (i.y)]  • 

The  last  statement  must  "be  qualified  in  two  respects,  First,  the  form 
(5)  applies  only  at  points  not  on  caustics  or  shadow  "boundaries.  Second, 
the  leading  term  will  "be  a  sum  of  expressions  of  the  type  given  in  (5) 
at  points  through  which  more  than  one  reflected  ray  passes.  There  will 
"be  one  such  summand  for  Roch  ray, 

♦  If  C  is  closed,  we  consider  the  solution  outside  C  only,  end  if 
C  extends  to  infinity  dividing  the  plane  into  two  parts,  we  con- 
sider the  solution  in  one  part  only. 
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If  we  apply  (2)  to  the  asymptotic  form  of  E  we  o"btain 

(6)  E'^'A^(x,y)  erp  [ilQ|;^(x,y)]  +  A^(x,y)  exp  [iTn|/^(x,y)]  -  0  on  C. 
IVom  (6)  we  have 

(7)  A^(x,y)  =  -  A^(i.y)     on  C 

(8)  M'j.(x»y)  =  Wj^i^fj)  on  C. 

(If  B  consists  of  a  sum  of  terms,  the  remaining  terms  occ\ar  in  pairs 

satisfying  (7)»(8),  where  the  indices  i  and  r  refer  to  respective 
terms  of  the  pair, ) 

A  and  >|r  off  C  can  be  determined  "by  the  methods  of  geometrical 

optics,  xising  the  values  on  C  given  by  (7),  (8),  In  particv-lar,  y 
se.tisfies  the  eiconal  equation 

(9)  1  vv^l  =  1. 

Once  \)f  is  found,  A  can  he  constructed  hy  the  conservation  of  energy 

formula  /  d  <r  I  ^^^ 

(10)  A^(x.y)  =  (  ^J  A^(x^,y^)  . 

Here  x,y  and  x   ,y  are  two  points  on  e  reflected  ray,  and  6.^  ^  d.  cr 

respectively  are  cross-sectional  lengths  of  an  infinitesimal  tube 
of  rays  at  each  of  these  points.  If  we  choose  for  x  ,y  a  point  on 
C,  then  (10 )  determines  A  off  C,  If  D  is  the  distance  along  the 
ray  from  x  ,y  to  x,y  and  if  K  is  the  curvature  of  the  reflected 

wavefront  at  x  ,y  ,  then  it  is  easily  shown  that 

On  the  basis  of  (5)-(ll)  the  asymptotic  form  of  the  re- 
flected field  may  be  immediately  written  as 

(12)   E^(x.y)'^  -(1  -  DK^)"'/^  \(^o»^o^  ®^  t^i^^o»^o^  "*■  ^^^   ' 

r 

The  point  x  ,y  is  the  point  on  C  from  which  the  ray  through  x,y  is 

reflected,  and  D  is  the  distance  between  these  two  points.  If  there 
is  no  reflected  ray  through  (x,y)  then  E  (x,y)  ***  0« 
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When  DK     =  1.  A    becomes   iirfinite.     This  occiirs  at  the  caustics  of  th« 
0  r 

reflected  rays. 


III  Ihe  Field  at  a  Caustic 

In  order  to  obtain  a  finite  value  for  E^  at  a  caustic,  we 
■begin,  with  Green's  formula 

(13)      E(x.7)  =  \(x.y)  +  i~- 

Here,D  is  the  distance  from  the  integration  point  x^ty^  oi  C  to  x.y^ 
and  H  '  is  the  zero-order  Hankel  function  of  the  first  kind.  To 

0 

obtain  the  asyniptotic  form  of  Fat  any  point,  and  in  particular  at 
or  near  a  catistic,  we  insert  into  the  integrand  the  asynrptotlc  form 

of  -P-,  which  can  he  computed  from  (4)  and  (12).   We  also  replace 

H^^^  "by  its  asymptotic  form  and  then  evaluate  the  integral  asymptotically. 

In  this  way  we  get  hack  the  expression  for  E^  given  "by  (12)  at  all 
points  off  the  caustic,  as  we  must  expect.  At  the  ca'ostlc,  however, 
we  obtain  a  finite  result  rather  than  the  infinite  value  given  by  (12). 

To  justify  the  method  described  we  must  assume  that  the 
caustic  does  not  intersect  C.  Porthermore  the  lower-order  terras  in 

the  expansion  of  P-,  which  have  been  neglected,  must  not  contribute 
terms  to  the  integral  which  are  of  the  flame  (or  higher)  order  as  the 
contribution  of  the  leading  terms.  That  they  do  not  can  easily  be 
proved  to  be  the  case  if  reasonable  hypotheses  are  made  about  the 
omitted  terms,  but  this  will  not  be  considered  here. 

We  will  now  carry  out  the  calculations  Just  described, 

V>,en  —  is  commuted  from  {k)   and  (12)  and  the  asymptotic  fom,  of 
Sn 

H^''''  is  employed,  equation  (13)  yields 

0 

(li^)   l!,.(x.y)  -  k^/2(Q„)-l/2  y  A^[^  -  ^-^/2exp[ik(v^4.D)-  ^]   ds . 
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35x18  integral  can  be  evaluated  asymptotically  "by  the  method  of  station- 
ary phase.  The  phase  \p.  +  D  is  a  fvmction  of  the  integration  point  on 

C,  and  can  he  considered  as  a  function  of  arc  length  s  along  C.  The 
points  of  stationsoy  phase  are  determined  "by 

.2 
'Then   if  — r-  (\)/  +  D)  ^  0  at  the  stationary  point,  (l4)  yields 
ds^ 


(16)   E^(x.y)  -  -^  (  ^  -  ^)  D-^/2  1  ^  (Vi  +  D)  1-^2 


exp 


ik^  +  D)  -  ^ 


1  -  sgn  — 2  (\|/^  +  VJ 

ds 


In  (l6)  all  quantities  on  the  right  side  are  to  "be  evaluated  at  the 
point  of  stationary  phase  determined  "by  (15).  If  there  is  no  station- 
ary point  the  right  side  is  to  "be  replaced  "by  zeroj  if  there  is  more 
thsn  one  stationary  point,  the  right  side  is  to  "be  a  sum  with  one  term 
of  the  a"bove  type  for  each  such  point. 

It  follows  from  (15)  that  at  the  stationary  points  either 
the  law  of  reflection  is  satisfied  or  the  reflected  ray  is  the  con- 
tinuation of  the  incident  ray.  In  the  latter  case,  evaluation  of 

(16)  shows  that  E  ~'  -«E.  so  that  E^x-'  0.   This  occurs  in  the  shadow, 
r    1 

In  the  former  case  (l6)  is  found  to  coincide  with  (12),  as  is  to  "be 

e3qi>ected. 
caustic. 


e3qi>ected.  It  is  to  "be  noted  that  the  phase  jumps  "by  «-  across  the 


^2 

The  eaqpression  (16)  is  not  valid  when  — ^  (t}r.  +  D)  =  0, 

But  if  for  a  point  x,y  this  second  derivative  is  zero  rnd  also  the 
first  derivative  (15)  is  zero,  then  it  can  he  shown  as  follov;s  that 
t'r.e   point  lies  on  a  caustic  of  the  reflected  ray,  \|/.  +  D  is  a 
fxiJictlon  of  T,y,?,  say  g(x,y,s). 
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For  each  value  of  a  the  equation  g„  (x.7.s)  =  0,  which  is  (15),  defines 
a  ray.  If  g   (xty»»)  =  0  also, .then  the  point  x,y  lies  on  the  envelope 

So 

of  the  family  of  rays,  which  is  the  caustic. 

In  order  to  determine  the  location  of  the  caustic,  we  compute 

.2 

S-—  (xfr  +  D)  and  obtain 

ds 

2 

(17)  i^  (x,,^  +  D)  =  (D"^  +  K^)  cos^oc   -  2Zg  coso(  , 
ds 

Here  K.  and  Z  are  the  CTirvatureo  of  the  incident  wavefront  and  of  the 
i     c 

curve  C,  respectively,  at  the  stationary  poiht.  The  angle  a     is  that 
between  ^.  and  the  normal  to  C  at  the  stationary  point.   The  express- 
ion in  (17)  vanishes  when  D  is  given  by  (see  also  j5j  ) 

♦  —1 

(18)  D  =  (2Z  -  K.  cos  o^  )    C08  0<  • 

Thus  (18)  determines  the  distance  D  from  C  along  a  ray  to  the  caustic, 
and  thus  the  curvature  Z^  of  the  reflected  wavefront  at  C  is  given  by 

V    -   (d*)"^  •   The  restdts  (l7J*0-8)  are  used  in  showing  that  (16) 

coincides  with  (12), 

To  obtain  the  field  on  the  caustic  we  must  evaluate  the 
integral  in  ilk)  differently.  For  this  purpose  we  consider  the 
integral  W  given  by 

(19)  W  =  y  A  (a)  exp  fifcf  (s)   ds . 

Stippose  that  at  a  stationary  point  «  at  which  f '    (s)  =  0  we  also 

have 
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(20)   f'(8)  =  f*'(7)«  ._   =f^^-^^  (7)  =  0, 


Then  the  Taylor  series  of  f(e)   becomes 


(21) 


(n). 


'(a)  «  f(T)  +  (a  -  F)^    hn^      * 


Therefore  W  la  asymptotically  given  by 

(22)    W  ^   A(r)  exp  [lkf(7)]     y  exp  [  Ik  ~7^ 
Thua  we  finally  obtain 


—Nil 


(8  -  a) 


(23)     W  ~  A(b)  e^  [  ikf  (B-)j  23c^  Pd  4  n'^)  in*)^  \   f  ^^^S") 


ds  , 


-1 

I  n 


where 


(24) 


F^=  exp  [(ijn)  8gnf^^)(8)j 


=  coa  "2  ^ 


if  n  is  even 


if  n  la  odd  . 


Let  us  apply  the  result  (23)  to  the  integral  (I'l)  when  the 
point  x,y  lies  on  a  caustic  and  when  the  first  n-1  derivatives  of 
tjr.  +  D  vanish  at  the  stationary  point.  Then  we  have 


(25)  E^  (z,y)'-T^2   n  (ni)n  |   (1  +  n~^)  (nD)^V^ 

•  exp  rik(M'i  +  D)  -  i  5]   7 

-'J 


ds 


-1 
a  coaor 


where 
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r  =  expl  (l5i)  Bgn  (v.  +  D)        if  n  is  eren 

n      f-       d  ^gii   1    J 

(26) 

P  =  COB  tt/2ii  if  n  is  odd. 


It  is  of  interest  to  notice  that  the  power  of  k  which  appears,  namely 

k^   ^  ,  is  always  a  positive  power  on  a  caustic  since  n  ^  3.  Ihus  the 
field  increases  indefinitely  as  k  does,  Ihe  larger  the  ralue  of  n  the 
higher  is  the  power  of  k,  as  one  should  expect  "because  as  n  increases  the 

phase  of  the  integrand  "becomes  more  constant.  If  the  phase  is  exactly 

1/2     V 

constant,  then  the  power  k  '   occuxIub,  as  will  be  seen  in  the  example  of 

the  para"bolic  cylinder, 

IV,  Comparison  With  the  De"bye  -  Picht  -  Lune'buX/g  Solution 

In  their  treatments  of  the  field  at  a  focus  or  caustic,  De"bye, 
Picht  Emd  Lune"burg  each  construct  an  exact  solution  of  the  wave  equation, 
or  of  Maxwell's  equations,  which  is  reguljrr  in  the  infinite  homogeneous 
space.  De"bye,  Picht  end  Lone  "burg  determine  their  solution  "by  requiring  it 
to  have  a  specified  "behavior  at  infinity,  and  Lur^eburg  proves  that  it  must 
reduce  to  a  specified  geometrical  optics  field  as  k  becomes  infinite. 

In  most  physical  problems  the  field  must  satisfy  boundary  condi- 
tions on  certain  surfaces,  and  satisfy  the  vmve  equation  or  Maxwell's 
equations  in  a  partly  inhomogeneous  medium.  These  conditions  are  not 
satisfied  by  the  Debye  -  Picht  -  Luneburg  solution.  However,  this  solution  does 
have  the  same  behavior  at  infinity  and  the  seine  geometrical  optics  field 
as  the  corresponding  physical  problem. 

The  question  therefore  arises  as  to  whether  the  Debye  -  Picht  - 
Luneburg  solution  implies  that  the  field  has  the  same  asymptotic  behavior 
at  a  caustic  as  that  of  the  corresponding  solution  of  the  physical  problem. 

By  the  considerations  of  Section  III  (first  and  second  paregraghs) 
both  of 'these  solutions, and  in  fact  any  solution  with  the  correct  geometri- 
cal optics  fieldjWill  imply  the  correct  behavior  at  caustics  and  foci  pro- 
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vided  that  the  lower-order  terms  in  the  integrand  of  (13)  cLo  not  contrilmte 
higher-order  terms.  Thus  all  Detye  -  Picht  -  lunehurg  solutions,  although 
represented  "by  less  convenient  integrals,  will  lead  to  the  same  leading 
asymtpotic  term  at  a  caustic  as  does  equation  (l4). 


V.  ApTilicatione;  1«  The  Faratolic  Cylinder 

Consider  a  cylinder,  the  cross  section  of  which  is  a  segment 
of  a  pejahola  (see  Figure  l).  Let  its  focus  he  the  origin  of  a  Cartesian 
coordinate  systeni,  its  axis  he  the  x-  axis,  and  its  equation  he 

(27)  y^  =  -  '+p(x  -  p)  x^a. 

As  is  Indicated  in  the  ahore  eq;'iation,  the  segment  of  the  parahola 
under  consideration  lies  to  the  right  of  x  =  a  and  is  concave  toward 
the  left,  since  we  assume  that  the  focal  length  p  is  positive. 

We  now  suppose  that  a  plane  wave  is  incident  from  the  left 
with  its  direction  of  propagation  parallel  to  the  x-  axis.  We  assume 
that  this  wave  is  given  hy 

(28)  E^  =  e^  . 

Tlie  incident  rays,  which  are  the  orthogonal  trajectories  of  the 
incident  wsvefronts  x  =  constant,  are  ohvioxisly  the  lines  y  =  con- 
stant parallel  to  the  x-  axis.  Of  these  rays  all  those  which  inter- 
sect the  paraholic  segment  will  he  reflected  through  the  focus,  and 
therefore  the  reflected  wavefronts  are  segments  of  circles  with  the 
focus  as  origin.  On  the  hasis  of  this  fact  and  equation  (l2),  we 
have  for  the  singly  reflected  field  E^ 
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Figure  1 
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(29)      B^^^(x.y)  -^  -  >/?.-*.!  e  ^^^^  "  ""^     ^^^''''^   ^°*^J 


E(^)(x.y)  -  ./IZl      e 


ik(2p  +  r)-lV2 


after  focus  , 


In  equation  (29)  r  ia  the  distance  of  the  point  1,7  from  the  focus 
and  D   is  its  distance  along  the  reflected  ray  from  the  point  of 
reflection  on  the  parabola.  The  phase  is  easily  determi-ied  since 
all  rays  reach  the  foCus  with  the  same  phase  ik2p» 

Some  rays  intersect  the  parabolic  segment  a  second  time 
after  having  heen  reflected  through  the  focus.  In  order  to  dis- 
tirtguiGh  these  rays,  \-re   introduce  the  angle  0  between  the  once- 
reflected  ray  and  the  z-  axis.  We  also  introduce  the  angle  y^ 
between  the  x-  axis  and  the  line  from  the  focus  to  tlie  iQjper  end  of 
the  peirabolie  segment.  If  y<  ^  then  a  >  0  and  the  entire  segment 
lies  t6  the  rieiit  of  the  focus.  In  this  case  no  ray  is  twice  re- 
flected. On  the  other  hand,  if  we  consider  the  complete  parabola, 
then  T  =  IT  and  every  ray  is  twice  reflected,  except  the  ray  in- 
cident along  the  x-  axis  and  characterized  by  ^  «  0,  In  all  other 

cases  (2*^  T  <  Ti)  all  rays  are  twice  reflected  except  those 
satisfying 

(30)      n-  T   >fi  >-(tt-  r ). 


which  are  once  reflected. 

After  the  second  reflection,  the  reflected  ray  tc-comeB 
parallel  to  the  x-  axis^ and  the  doubly  reflected  field  E^  ^  (x,y) 

is  found  from  the  equations  (12)  and  (29)  to  be  given  by 


(31)    z^l\^.7)    -  |£   .  "=^'*  -  ^'-'"'^ 


y 


l^ 
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The  cojnplete  geometrical  optics  field  can  now  te  descrited 
as  zero  in  the  shadow,  E.  at  points  outside  the  shadow  throti^  which 

no  reflected  ray  passes,  E,  +  E^  'at  points  throtigh  which  one  re- 

1    r 

fleeted  ray  passes,  B  +  E^  ^  +  Ej;i^  at  points  through  which  two 

singlj'  reflected  rays  pass,  E,  +  E^  '  +  E^  'at  points  through 

which  one  singly  and  one  douhly  reflected  ray  pass,  and 

\  *  r   "•■  ^  «  "♦■  ^r  ^^   points  through  which  two  singly  and  one 

doubly  reflected  ray  pass. 

We  now  observe  that  the  field  E   ,  as  given  by  (29), 

becomes  infinite  at  the  focus  of  the  parabola,  and  that  the  field 

(2) 

Ejl   as  given  by  (31)  becomes  -infinite  along  the  axis  y  =  0.  This 

latter  singularity  appears  only  in  the  case  of  the  conrplete  parabola, 
since  the  doubly  reflected  ray  y  =  constant  s\tffers  it  first  reflec- 

ion  8/  a  point  on  the  parabola  with  the  ordinate  — ~  .  Thus  as  y 

y 
approaches  zero,  the  ordinate  of  the  point  of  first  reflection  nuat 

become  infinite,  and  this  is  possible  only  with  the  cosplete  parabola  • 

The  X-  ezis  is  not  a  caustic  in  the  ordinary  sense  of  an  envelope  of 

rays,  but  the  geometrical  optics  field  becomes  infinite  there  becaiose 

the  energy  between  incident  rays  with  an  arbitrarily  large  separation 

is  reflected  between  rays  with  an  arbitrarily  small  separation. 

We  will  now  assume  that  the  parabolic  segment  is  finite 

(V  <  tt),  and  we  will  obtain  a  finite  expression  for  the  field  at 

and  near  the  focus,  which  is  the  only  caustic  occiirring  in  this  case* 

To  this  end  we  employ  (1^),  but  we  need  not  consider  the  contributions 

(2) 

of  E    to  the  integrand  since  the  doubly  reflected  field  contains 

r 

no  ray  through  the  focus  and  therefore  contribtites  only  loweivorder 
terms  to  the  integral. 

First  we  introduce  the  parabolic  coordinates  f  »  ^  defined  by 

^2  =  (x2*y2)2  .  ,.      ;^  2  =  (x^  ^  yV  ^, 


(32) 
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In  these  coordinates  the  eq-uation  of  the  parabolic  segment  is 


(33)        yi  =  tV^  .  0  ^  ^   ^  y  2(p  -o()  . 


Nov   (l4)  hecomes 


i3k)     E^^^(i.l?)~-/^    espdkp-l 


4 


V2(p-o<) 


r-1 


y 


2(p  -«<) 


-1 
2 


D  2  esp(i)d)-    I  $   2^ 


^f 


=  -y2p 


At  the  focus  this  yields 


(35)   E^^^0.0)'^'-2)^  e:cp(2ilq)-ij) 


V2(p-«') 


(Ho'^p')'  *f 


^.2/^.zp(2iki>-if) log  /i-^  +/775:  I  . 


To  obtain  the  field  near  the  focus  we  introduce  the  polar 

coordinates  P  ,  Q  with  the  focus  as  origin  and  with  p  defined  hy 
^  1 

P  =  k(x2  +  7  )2  >  Thus  p  is  2n  times  the  distaiace  measured  In 
wavelengths,  of  a  point  from  the  focus.  Now  (3^)  becomes 


-Ik- 


y 


(36)         E^^\p.«)  _  ./S    e.5>(2  ikp  -  ^)^^  sec  |  •  e^P  ^^^^^^^  d/5 


Here  V  =  cos  ~     j]  2!^^  3   •'•^  *^®  angiilar  coordinate  of  the  end  of 


2lp— ct, 


the  paraholic  segment,  ffigures  3-^»  "based  on  eq,  (36),  show  the 
variation  of  the  reflected  field  eiqDlitude  |E^  |  along  and  at  rj 
angles  to  the  axis  for  two  different  parabolic  segments. 


2,  The  Circular  Cylinder 

A  plane  wave  exp(ilac)  is  incident  on  the  concave  side  of  the 
sector  of  a  circular  cylindrical  mirror  of  radius  R.  The  situation  is 
represented  in  Hgure  5»  Ihe  center  of  the  mirror  has  heen  plaxsed  at 
the  origin  of  a  cartesian  coordinate  system  so  that  the  incident  wave 
moves  parallel  to  the  z-a:!:is  in  a  positive  direction.  The  angular 
extent  of  the  mirror  is  a  +  3,  where  a  measxires  the  angle  "between 
the  z-asis  and  the  end  of  the  sector  in  the  first  quadrant,  and  3 
measures  the  angle/  "between  the  x-asis  and  the  end  of  the  segment 
in  the  fourth  quadrant. 

We  shall  consider  only  sectors  such  that  no  ray  is  reflected 
more  than  once.  It  is  apparent  from  Figure  5  that  P  can  he  at  most 
TT  -  3a, 

The  phase  of  the  incident  wave  at  a  point  of  reflection 
on  the  mirror  is  R  cos  ©,  where  6  is  the  angle  of  incidence;  from 
Figure^ it  is  apparent  that  the  equation  of  the  reflected  ray  is 

(37)    (R  sin  0  -  y)  /  (R  cos  e--  r)  =  tan  20  j 
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Figure  5 
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(x,y)  is  an  arbitray  point  on  the  reflected  ray. 

Let  D  "be   the  distance  the  reflected  va^efront  has  traveled 
from  the  point  of  reflection  meastured  along  the  reflected  ray.  Dien 
the  phase  ^   of  the  reflected  wavefront  at  the  point  of  ohservation 
la  D  +  B  cos  Q,     It  is  clear  from  Fig.  5  that 

(38)  D  =  (E  cos  «  -z)  /  cos  20  . 

Thta  the  total  phase  at  the  point  of  ohservation  is  given  "by 

(39)  \y  =  E  cos  «  +  (E  cos  Q  -  z.)  /   cos  20  • 

By  a  similar  consideration,  ©r  from  (37 )»  we  obtain  an  expres- 
sion for  \f   involving  y.  Using  this  eacpression  and  (39)  we  obtain  the 
following  parametric  equations  for  a  reflected  wavefront  i.e.,  a  curve 
of  constant  phase  v)/  : 


X  =  (B  cos  0  -  \p)  cos  20  •♦•  E  cos  0 

(^0) 

y  =  (E  cos  0  -  \jr)  sin  20  +  E  sin  0  . 


We  can  tise  (^O)  to  obtain  a  coordinate  system  for  points 
(x,y)  on  a  reflected  ray  in  terms  of  the  distance  D  and  the  angle  of 
incidence  0  by  replacing  ly  by  E  cos  «  +  D: 


X  =  -  D  cos  20  +  E  cos  0 

(^1) 

y  =  -  D  sin  20  +  E  sin  0  • 
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Now  we  ©"bserve  that  the  Jacobian  of  this  trajQsf  onviation  is 
2D  -  H  cos  Q,   which  is  zero  when  D  =  ~  R  cos  «.  This  is  the  distance 
D  to  the  caustic  given  "by  our  formula  (18): 


(42)    D*  =  (2k:^  -  k^  cos  Q)'^   cos  Q  =  (Z/'R)'^   cos  «  =  | 


B  cos  4  • 


We  can  check  this  result  hy  calculating  the  curvature  of  the  reflected 
wavefront  vjr  =  constant  froa  (4o);  we  obtain  the  result  Z=   2(E  cos  9)~ 
as  expected,  Fron   our  fornialas  (12)  we  have  for  the  reflected  geometrical 
optics  field  the  q\iantity 

(43)    E^ (1  -  5-3^-5  r^^^    exp  [  ik(E  cos  «  +  D)  ] 

in  terms  of  the  coordinate  system  defined  "by  (4l), 

At  the  caustic,  where  D  =  t'  R  cos  9,   the  reflected  field 
takes  a  different  form,  which  is  found  from  (25): 

1-1  i         3  l~- 

ikh)       ^^-^  ik3)l -n 'Z    P  (i)  cos  f  .  k^  v^cosT-  ^^  (Eco8«  +  D)  |  \xp[|ikacos«-i^  , 

dO 


To  calculate  the  third  derivative  in  this  espression  we  make  use 
of  the  fact  that  ^  =  0  at  the  point  of  reflection  and  that  D  =  ^  R  cos  « 

there.  Then 

(45)  r^  8in(0  -  «^)  =  J  R  sin  2«  . 

In  addition  we  use  the  fact  that  — ^  =  0  when  the  observation  noint  is  on 
the  caustlCj  and  obtain 

(46)  r^  cos(«  -  «^,)  =  I  (1  +  sin^  «)  , 
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Then,   -"asing  (45)  and  {^6)^   it  is  easy  to  find  that  when  the  point  of 
observation  is  on  the  catistic, 

(47)  ^   =  -  6E  sin  «  . 
dO-^ 

Hlhtis  the  reflected  field  at  the  caustic  Is 

-1  11-1 

(48)  E^  '^  2n^  P  (|)  cos  g  •  (35E)°  cos^  «  .  sin^  «  ezp[  |  ikR  cos  «  -  i  J]  , 


TiiiB   result  holds  everywhere  on  the  caustic  except  at  *  =  0,  where  the 
caustic  has  a  cusp.  To  obtain  the  field  at  this  point  ve  need  one  more 
derivative  of  \jr  s.t  the  cusp;  this  is  given  by 

(49)       ^   =  -  6E  . 

Bie  reflected  field  at  the  cusp  is  then 
1-1  i 


(50)  E^-'  2^  v^  r  ih   •  (kR)^  exp  [  I  ikR  -  ^  iTt]   . 


In  the  neit^hborhood  of  the  caustic  curve  there  is  transition 
^   between  the  various  forms  of  the  reflected  field  given  by  (44),  (48)  aiadC^t)) 


.  lo  obtain  this  transition  form  we  must  return  to  the  integral 
expression  from  which  the  general  results  were  derived.  In  fact  we 
write  for  the  reflected  field 

-1  1-1 
(51)  Ej.^  (2rr)^  k^  D^  E  cos  «  exp  [  ik(H  cos  «  +  D)  -  i  Jj  /"exp  J  d« 


where 
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(52)  $  =  |ik  I  (D-^  R^cos^*  -2R  co8«)  (««-«)^+  ^"^  Hasina*  (l^-EcosO)  («•-  O)^ 


£:xcept  in  the  neighborliood  of  the  cusp.  In  the  neishhorhood  of  the 
cusp  we  have 


2.    r  1  -^-2  T)2_4_oA   /T\_T3/>«ori^1    ^n«_  fl^ 


|)=:  ik  J  [|d~^  E^cos^*  -Rcos«]    («•-  «A   [5D-2  E23i^2«  (D-EcosO)]    («•-  O) 


(33J  +  |.  [eco3«  -D"^  R  (B^Dcoad)  +  klT^  R^sin^  «  -3D"^  H^  (R-Dco8«)' 

+  18  D"^  R^8in^«   (H-Dco8«)  -  ISD*"^  R^8in^«  ]    («'-«)        j-    . 


If  we  vse   the  value  of  ^   given  hy  (53) »  the  integral  in  (51 )  can  he 
evaluatec  easily  in  terms  of  standard  functions  in  the  special  case 
where  the  ohservation  point  i8  the  ray  reflected  throu<f^h  the  cusp,  the 
ray  heing  given  "by  9=0.  On  the  other  hand,  if  we  use  the  value  of 
^        given  "by  (52),  we  can  evaluate  the  integral  in  (51)  with  no 
restriction. 

Consider  the  integral  in  (51)  witL  §    given  hy  (52).  It  has 
the  form 


.»  -1 


(54)   I»  /  eip  [i  (rt3+d^t2)3  dt-  lr|  ^/  e^  [i  (8^  +  /r  ^  g^)]  d8. 


^V-2/3 


Let  8  =  u  -  * — r .  Then 
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—1                                                       3         2 
I  =  2|Tpexp   [^]/     cos    (  lJ73-)du. 

(55)  -13  2 

=  2|rPe:cp  r^l    Ci,  ( ^^7^  )     . 


where  Ci  (a)  is  Hardy's  generalized  Airy  integral  given  "by  Watson 
n 


3         3 
According  to  (52)  we  have 


Y"  =  (kE^  sin  2©  /  ^D^)  (D  -  E  cos  Q) 

(57) 

/  =  I  k  [  R^cos^  e  /  D)  -  2E  cos  «  ]   , 


and  for  the  reflected  wave  field  in  the  neighborhood  of  the  caustic 

(58)  E^'>/(2k/Tffi)^|  rl^  cos  Q   exi.rik(E  cos  «  +  D)  +  ±{^=^—^^  r-)3Ci,(  '^/^   ) 

It  is  understood  of  course  that  (58)  holds  only  when  the  point  of  oh- 
servation  is  away  from  the  cusp. 

If  the  point  of  observation  is  in  the  neighborhood  of  the  cusp 
and  on  the  line  *  «  0,  we  can  use  (51)  and  (53)  to  obtain  the  reflected 
wave  field  in  a  form  similar  to  (58).  Here  the  integral  on  the  right 
of  (51)  has  the  form 
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(59)  I  =     /     ezpf  -i(  coq'V  P  «'2)    I   d«'  =  a  ^/  exp   [-i(tV  Pd  ^  t^)  ]   dt 

1/-0D  '/-oo 

We  assTime  a>  0  in  the  region  considered,  since  D'^^  there.     Then 

-1  p  -1  -1 

(60)  I  =  aa  ^  exp   (iP^Sa)  J  C^^^    (^  Pa  2)  -i  S^^  (J  f^^a  2)   I  , 


Thus 

1-1  (  -i  -1    "j 

(61)  Ej.-'(2k/nD)2  2c,  ^  exp   [ik(E  +  D)  +  i(P^/8a  -  tt/^)J      C^^^^J^^  ^)--^HU^^  ^)   [ 

with 

(62)  CL  »  If  [  BD~^   (E  -  D)  -  E  +  3R^  d"^  (E  -D)^  ] 

12-"' 
p  =  E  -  I  E^  D  -, 

igF.in  the  fiinctions  C  .    and  S  .    sxe  Hardy's  generalized  Airy  integrals'-  -^    : 

1  r 

C^^(a)  =  m^   [2n  sin   (Tr/Zn)]"-     J   J_^/^   (20^^)  -J^/^   (2a^/2) 

(63)  , 
S^^(a)  =  m^   [2n  cos    (tt/2q)]"^     .    J_^/^   (2a''/^)  +J^/^   (2a^/^)    I     . 

Eqiiation  (61 )  gives  the  reflected  field  in  the  neighlsorhood  of  the  cusp 
provided  that  the  o'bservation  point  is  on  the  axis  of  the  circular  mirror, 

AcknowledgEients 

The  authors  wish  to  thank  Professor  Fritz  Eeiche  for  his  stimulat- 
ing discussions  duringpiis  investigaticn. 


-21- 


Eeferences 


1.  Iiune"burg,  E.Z. ,  Asymptotic  Development  of  Steady  State  Electromagnetic 
Fields 

N.T.U.,  Washington  Square  College,  Mathematics  ReseeTch  Group, 
Besearch  Report  No.  EM-ll*.  (19^9). 

2.  Detye,  P.,  Das  verhalten  von  Lichtwellen  in  der  nShe  eines  Brennpunktes 
Oder  ciner  Brennlinie 

Ann.  D.  Phys.  (4)  21,   755,  (1909). 

3.  Picht,  J.,  Tiber  den  Schwingungsvorgang,  der  einem  heliehigen  lastigmati- 
schen)  Strahlenlnlndel  entspricht 

Ann.  d.  Phys.  (4)  22.  685,  (1925). 

ii,   luneturg,  R.Z. ,  Asymptotic  Evaluation  of  Diffraxition  Integrals 

N.Y.U, ,  Washington  Square  College,  Math.  Res.  Group,  Research 
Report  No.  SM-15, (19^9). 

5.  Keller,  J.B.  and  Keller,  H.B.,  Determination  of  Reflected  and  Transmitted 
Fields  hy  Geometrical  Optics 

N.Y.U. ,  EM;-13,  (19^9)  or  Jour,  Optical  Soc.  Am.,  Vol.  ijO,  No.  1, 
pp.  48-52,  January,  (19^9). 

6.  Watson,  G.N,,  A  Treatise  on  the  Theory  of  Bessel  Functions 

The  Macmillan  Co.,  New  York,  (1948),  p.  320  ff. 


Date  Due 

lo  l^^hi 

.t^ 

wTSr' 

''-■'V  '  -^ 

1   "^ 

■&f^3   m 

DEC  1  5  1910 

MAY  10  • 

08S 

f) 

PRINTEC 

IN  U.  S.  A. 

of  toe  field  at  a 


sM-55 


Kay 

Asymptotic  evaluation  of 

the   field  at  a  caustic. 


I 


